Orbital selective insulator-metal transition in V2O3 under external pressure 
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We present a detailed account of the physics of Vanadium sesquioxide (V2O3), a benchmark system 
for studying correlation induced metal-insulator transition(s). Based on a detailed perusal of a 
wide range of experimental data, we stress the importance of multi-orbital Coulomb interactions in 
concert with first-principles LDA bandstructure for a consistent understanding of the PI-PM MIT 
under pressure. Using LDA+DMFT, we show how the MIT is of the orbital selective type, driven 
by large changes in dynamical spectral weight in response to small changes in trigonal field splitting 
under pressure. Very good quantitative agreement with (i) the switch of orbital occupation and (ii) 
S — 1 at each V 3+ site across the MIT, and (Hi) carrier effective mass in the PM phase, is obtained. 
Finally, using the LDA+DMFT solution, we have estimated screening induced renormalisation of 
the local, multi-orbital Coulomb interactions. Computation of the one-particle spectral function 
using these screened values is shown to be in excellent quantitative agreement with very recent 
experimental (PES and XAS) results. These findings provide strong support for an orbital-selective 
Mott transition in paramagnetic V2O3. 



PACS numbers: 71.28+d,71.30+h,72.10-d 
I. INTRODUCTION 

Correlation driven metal insulator transitions have re- 
mained unsolved problems of solid state theory of elec- 
trons in solids for more than five decades. The pioneer- 
ing work of Mott, 1 and of Gutzwiller, 2 Kanamori, 3 and 
Hubbard 4 involved a detailed exposition of the view sug- 
gesting that description of such phenomena lay outside 
the framework of band theory. Subsequent, more recent 
developments, including discovery of high-T c supercon- 
ductors, rare-earth based systems, colossal magnetore- 
sistive oxides along with whole families of other systems 
have shown that strong electronic correlations give rise to 
widely unanticipated, fundamentally new types of metal- 
lic behaviors, namely, non- Fermi liquid metals. A host of 
very careful studies now clearly show that these anoma- 
lous responses seem to be correlated with the existence of 
a correlated metallic state on the border of a Mott insula- 
tor in e?-band oxides, 5 or to a metallic state in proximity 
to a localisation-delocalisation transition of /-electrons 
in rare-earth compounds. 6 

The corundum lattice-based transition metal oxide sys- 
tem vanadium sesquioxide (V2O3) has been of interest for 
more than five decades as a classic (and by now a text- 
book) case of an electronic system with S — 1/2 local 
moments at each site exhibiting the phenomenon of the 
correlation-driven Mott-Hubbard metal-insulator transi- 
tion. 1,4 Widely accepted wisdom has it that this is one of 
the few cases where modelling in terms of a simple one- 
band Hubbard model is appropriate. Over the last few 
years, experimental and theoretical work has forced a re- 
vision of this view, leading on the one hand to a spurt of 
new approaches, and on the other to an expansion of our 
perception of what is new and important in the physics of 
TM-oxides in general. Specifically, taken by themselves 
as well as in combination with properties of other sys- 



tems like ruthenates, CMR manganites, etc, these new 
studies force one to refocus the attention in terms of the 
strong coupling and interplay of spin and orbital degrees 
of freedom (not to be confused with the usual spin-orbit 
coupling, though that may also be relevant in some situ- 
ations) and of their combined influence on the nature of 
charge and spin dynamics in TM oxides. 

In what follows, we aim to present a far from complete 
view of the questions which are posed by these new stud- 
ies. Focusing our attention to the early vanadium oxides, 
V2O3, we will start with a somewhat detailed perusal of 
earlier work, review key recent experiments, and follow 
them up with a discussion of their implications for the- 
ory. Finally, we will suggest a rather detailed scenario for 
correlation-induced metal insulator transitions in V2O3 
that ties together essential experimental constraints in 
one picture. In doing so, we will give a detailed descrip- 
tion of our theoretical modelling using a combination of 
the local density approximation (LDA) combined with 
multi-orbital dynamical mean field theory (DMFT) using 
the iterated perturbation theory (IPT) as the "impurity" 
solver. 

In the first part, we will confine ourselves to summa- 
rizing known -and not so well-known experimental re- 
sults on the effect of external pressure and Cr doping 
on the thermodynamic and transport properties of V2O3 
along with the magnetic and orbital structure (and their 
changes) across the metal insulator transition. 

In the second part, we will first review the earlier the- 
ories for the MIT in terms of the one-band Hubbard 
model as well as the more recent multiband Hubbard 
model (which allows a description in terms of a S = 1 
model). Finally, we will propose a new scenario: one 
where the abrupt change in the character of spin (and 
presumably also orbital) correlations across the MIT is 
described within the strong correlation scenario. In par- 
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ticular, we will show how a two-fluid (i.e., orbitally se- 
lective) description can be derived from first principles, 
and demonstrate how the properties of V2O3 can be un- 
derstood in this scenario. 



II. EXPERIMENTAL REVIEW 

The magnetic structure of V2O3 has been measured a 
long time ago. 7 Its interpretation has however remained 
a subject of controversy. The robust aspects are: the 
antiferro-insulator (AFI) is characterized by AF order 
which spontaneously breaks the crystal symmetry of the 
corundum lattice -and in modern parlance, it corresponds 
to the C-type AF order 8 with one F and two AF bonds 
in the hexagonal plane. The vertical V-V pairs (with 
V 3+ ) form dimers in the solid; these are aligned anti- 
ferromagnetically with the inplane V-V pairs, and fer- 
romagnetically with other vertical V-V pairs. In terms 
of these dimers, the corundum lattice can be viewed as 
a distorted simple cubic lattice. Based on the S = 1/2 
picture of Castcllani et al., the spin waves were "charac- 
terized" in terms of a Heisenberg-like model. 8 This was a 
commonly accepted picture for almost two decades, un- 
til recent experimental results forced one to reanalyze it. 
These are: 

(i) Change of magnetic correlations across the AFI to 
AF-metal (AFM) phase transition in V2- y C>3. In a one- 
band Hubbard model scenario, one would expect that 
the magnetic correlations in the AFM should be a rem- 
nant of those in the AFI, for c.g, broadened spin-waves of 
the AFI. However, measurements revealed, surprisingly, 7 
that the transport and thermodynamics is due to onset 
of magnetic order which seems to be totally unrelated to 
that in the AFI. The AFM is characterized by incommen- 
surate order with Q||c, in contrast to that for the AFI, 
characterized by Q = (1/2, 1/2, 0) (hexagonal notation). 
This order in the metal is reminiscent of a small-moment 
SDW derived from a Fermi surface instability. 7 Constant 
energy scans in INS also show that the magnetic fluc- 
tuations are far from conventional spin-waves; they are 
more reminiscent of particle-hole modes with extremely 
small correlation length of about 14 A, much smaller than 
that characterizing the spin- waves in the AFI. Finally, 
the AFI-PI transition (P meaning paramagnetic) occur- 
ing around T c = 150K is second order, once again accom- 
panied by an abrupt change in short-ranged spin corre- 
lations. An interesting correlation concerns the change 
of lattice structure from monoclinic (AFI) to corundum 
(PI,PM), and may be linked to a change in orbital cor- 
relations. 

Further, abrupt jump of the crystal volume (without 
change of symmetry) is also observed across the PI-PM 
transition at higher T. 9 The c-axis distance decreases 
abruptly at the PI-PM transition, while the a-axis dis- 
tance slightly increases, and this change needs to be cor- 
related with the conductivity jump at the transition (see 



below) . 

— > This abrupt switching of the spin correlations across 
the AFI-AFM transition is inexplicable within an ef- 
fective correlated one-band scenario. The AF-P transi- 
tion^) cannot be considered as order-disorder transitions 
of the usual type. 

— > Are orbital degrees of freedom involved in the mon- 
oclinic to corundum structural change at the MIT? 

— > What is the specific relation between the lattice 
contraction and the jump in the dc conductivity at the 
MIT? Does it involve changes in carrier concentration, or 
in the carrier mobility? 

(ii) Recent X-ray experiments by Park et al. have 
revealed the existence of an admixture of (e^e^) and 
(eg i7 aig) with i = 1,2 and a spin S = 1 on each V 
site, in contrast to the 5 = 1/2 proposed in Ref. 8. In 
addition, large differences in their ratio have been found 
in the AFI, PI and PM phases. In particular, this ratio 
is (e^,e^ 2 ) : (e^,a lg ) is 2 : 1 in the AFI, 1.5 : 1 in the 
PI and 1 : 1 in the PM phase. 

The AFI phase is also characterized by a monoclinic 
distortion involving a uniform rotation of all the V-V 
pairs, 11 an observation which puts specific constraints on 
possible orbital order in the AFI. However, the corundum 
structure is recovered in the PI phase. (Does this change 
involve switching of the orbital correlations across AFI- 
PI boundary, and if so, how?) Based on this, Ezhov et 
al. 12 have proposed a S = 1 model without orbital de- 
generacy as an alternative starting point to Ref. 8. More 
studies along this line have been done by Tanaka. 11 An 
alternative point of view by Shiina et al. 13 pictures the 
AFI as a C-type AF ordered state with S = 2 and ferro- 
orbital ordering. 

— > Starting point should involve S = 1, suggesting use 
of multiband models is necessary to describe V2O3. In 
orbitally degenerate cases, the ground state is simulta- 
neously spin and orbital ordered, and the strong cou- 
pling between the elementary excitations involving both 
spin and orbital flips is expected to result in emergence 
of qualitatively new behavior (not, however, if the sce- 
nario of Ref. 11 is taken to be valid). It is important to 
notice that even in the para-orbital/magnetic state(s), 
additional strong scattering resulting from coupled spin 
and orbital fluctuations might cause pronounced devia- 
tions from expectations based on (correlated) Fermi liq- 
uid theory. 

— > Notice that this change of orbital occupation also 
implies an important role for the trigonal distortion (this 
would act like an external field in orbital pseudospin 
space) . This quantity determines the occupation of the 
relevant orbitals, and in a coupled spin-orbital system, 
determines the effective exchange interactions (and hence 
the magnetic structure). In fact, the importance of 
this quantity, and its pressure (strain) dependence, has 
been identified in many members of the corundum based 
oxides 14 (FC2O3 in connection with the Morin or the 
strain-induced spin-flop transition, Cr203, in the same 
connection, Ti 2 03 as manifested by an anomaly in the 
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T-dependence of the A\ g phonon frequency across the 
MIT), suggesting a related common origin of the varied 
manifestations observed in this structurally related class 
of oxides. 

— > It might be interesting to look at possible anomalies 
of the Aig-mode frequency in Raman scattering measure- 
ments across the P-MIT in V2O3. 

(Hi) Resonant X-ray scattering (RXS) measurements 
at the V K-edge have been performed by Paolasini et al. 15 
The resonant Bragg peaks observed in the AFI phase are 
interpreted in terms of a 3d orbital ordering. However, 
the authors of Ref. 15 interpreted their results within the 
picture of Castellani et al, leading to difficulties with (i) 
and (ii). 

— > One requires a re-interpretation in a way consistent 
with (i) and (ii), for e.g, see Mila et al. An Aside: Is the 
Kugel-Khomskii type of modelling required for the AFI? 

(iv) Finally, Lovesey et al. 16 argue that the resonant 
Bragg peaks arise due to ordered orbital magnetic mo- 
ments of the V ions. Indeed, a large orbital contribu- 
tion Ml I Ms — —0.3 to the total magnetic moment was 
claimed in Ref. 16. If correct, the effect of spin-orbit cou- 
pling might become important, as argued in Ref. 11. On 
the other hand, if the spin-orbit coupling leads to small 
effects on the electronic structure, a re-interpretation 
would be in order. 

— > How important is the role of spin-orbit coupling in 
the AFI phase? 

(v) Optical spectroscopy provides a detailed picture 
of charge dynamics. Careful studies by Thomas et al. 17 
reveal all the characteristics of a strongly correlated sys- 
tem: an "upper Hubbard band" (UHB) feature with a 
threshold in the AFI (and PI), and a sharp, quasicoher- 
ent feature, along with an intense mid-infra-red peak and 
the remnant of the UHB on the PM side. Calculations 
within the framework of a S — 1/2, one-band Hubbard 
model 18 claim to obtain very good agreement with the 
observed spectra. However, to achieve this, the Hubbard 
U has to be changed by a factor of 2 on going from the PI 
to the PM state, which is hardly conceivable. Alterna- 
tively, the effective hopping, or the degree of itinerance, 
should increase in the metallic state, leading to a change 
in the effective U/t value, and, as is ubiquitous in strongly 
correlated scenarios, to a large transfer of spectral weight. 
In particular, given (i) — (Hi), this should involve carriers 
coupled to spin-orbital degrees of freedom along with the 
concomitant lattice distortion. 

— > What is the specific nature of the correlation be- 
tween the change in spin-orbital correlations across the 
MIT and the tendency to increased itinerance which 
drives this transition? In particular, is the abrupt change 
in the e^e^ : e? gi a\ g ratio related to increase in carrier 
concentration or the carrier mobility (kinetic energy)? 
Notice that (see below) the dc transport data can be rec- 
onciled with increase in the carrier density, 20 so this is 
an important point deserving more attention. Obviously, 
to make a plausible correlation between the two requires 
experimental characterization, as well as a proper treat- 



ment of these coupled correlations. 

(vi) Photoemission Spectroscopy (PES) 21 reveals fur- 
ther proof of the correlation-driven character of the 
insulator-metal transition. At high T, the data is claimed 
to be consistent with a "thermally smeared quasicoher- 
ent" peak, something within reach of single-site theories. 
However, it is more conceivable that strong inelastic scat- 
tering from coupled spin-orbital excitations gives rise to 
non-quasiparticle dynamics in the PM phase. At lower 
T, appreciable changes are observed in the PES spec- 
tra across the AFI/PM and PI/PM transitions, with 
the characteristic transfer of spectral weight from high 
to low energy over a scale of almost 4 eV (notice that 
Tmi — 300K) implying a drastic rearrangement of elec- 
tronic states over a wide energy scale. 22 In the PM state, 
the PES spectrum shows the asymmetric two-peak struc- 
ture with a Fermi edge (but we draw attention to the fact 
that this low-energy "peak" is anomalously broad, sug- 
gesting non-quasiparticle dynamics), while a clear open- 
ing of a spectral gap (E g ) occurs in the AFI and PI 
phases. In the AFI phase, Ef FI ~ 0.3eV, while for the 
PI, E^ 1 ~ 0.23eV, with a more symmetric lineshape. 

Earlier PES studies across the PI-PM transition have 
been controversial; in particular, the question of the In- 
dependence of the low energy spectral weight was not set- 
tled till recently. Quite recently, this question has been 
answered by the Michigan group, 23 and the T-depcndcnt 
renormalized, and heavily damped "quasiparticle" con- 
tribution has indeed been observed. 

Details of the PES spectra at high T are seemingly well 
captured by a DMFT applied to a multiband Hubbard 
model in combination with the actual LDA bandstruc- 
ture (see below for more details). 22 There are still some 
discrepancies between LDA+DMFT and experiment at 
lower T < 400 K, however: the "quasiparticle peak" is 
too broad by a factor of 2 — 3, and the details of the 
PES lineshape in the PI still remain to be calculated. 
It is possible that screening-induced renormalization of 
U, etc. needs to be included; however, it is a very dif- 
ficult task to do this from an ab initio starting point. 
Alternatively, or in concert with the above, the dynam- 
ical effect of intersite correlations might be expected to 
become increasingly important at lower T. Such effects 
are out of scope of LDA+DMFT, and require extensions 
to treat dynamical effects of spatial correlations, a more 
demanding task. 

Given this, the interpretation of the PES spectrum 
will also need a rc-cxamination. In particular, the possi- 
ble importance of short-ranged spin-orbital correlations 
might be required to understand the anomalously broad, 
low-energy feature observed in PES in the PM phase. 

In a coupled spin-orbital system, the degree of itiner- 
ance is directly related to the changes in spin and orbital 
correlations coupled to possible structural changes. In 
course of its hopping motion, an electron is scattered by 
coupled spin-orbital excitations, i.e, by simultaneous flip- 
ping of spin and orbital pseudospins. In the AFI, this is 
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not enough to destroy AF/O order (O meaning orbital). 
An understanding of the change in AF/O correlations 
across the MIT is necessary to understand the enhance- 
ment of itinerance. In particular, within the framework 
of Shiina et at, do the changes involve a cooperative melt- 
ing of the AF/O order of the AFI? How does one then 
try to understand the PI/PM transition? 

(vii) One of the most spectacular hallmarks of the I-M 
transition in V2O3 is the sharp jump in conductivity by 
seven orders of magnitude! Is the jump of a(T) 19 driven 
by a jump in the carrier density at the transition, or by 
an increase in the mobility? 20 Hall effect measurements 
would be a probe to answer this question (complications 
due to possible relevance of spin-orbit coupling, if im- 
portant). On the barely metallic (close to the AFI) side 
(V2- y 03), the Hall constant Rh{T) shows behavior rem- 
iniscent of the cuprates; 24 it is strongly T-dependent, in- 
creasing with decreasing T with a peak around the AF 
ordering temperature, followed by a drop at lower T. 
The T-dependence gets weaker with increasing metallic- 
ity (y)- More similarities with the normal state of the 
high-T c cuprates are seen in the different T dependences 
of p(T) ~ T 3 / 2 and cot6 H (T) ~ aT 2 + b for small y, 
which evolves into more conventional FL behavior with 
increasing y. 24 Such a behavior would mandate strong 
local moment scattering in the metallic phase. Given the 
strong correlation signatures observed globally, a descrip- 
tion in terms of vagaries of the Fermi surface is untenable. 

— > Docs the I-M transition involve a jump in the carrier 
density? 

— > How does one understand the anomalous features 
of the Hall data near the MIT? In particular, similar- 
ity to cuprates suggests that such anomalies might be 
more general manifestations of the breakdown of Fermi- 
liquid theory (FLT) near the Mott transition to a Mott- 
Hubbard antiferromagnet, as opposed to a Slater anti- 
ferromagnet. The observation of overdamped spin waves 
with extremely short correlation length and anomalously 
broad linewidth is also reconcilable in terms of a strong 
scattering scenario. 

A. Summary of experimental results 

In conclusion, experimental results reveal very inter- 
esting points concerning the nature of the ground states 
and collective excitations in the different phases of V2O3. 

1. In the AFI phase 

(i) C-type AF order with ferro-type concomitant or- 
bital order. In terms of the V-V pairs, it corresponds 
to C-type AF order on a distorted simple cubic lattice. 
A Kugel-Khomskii type of model is required to derive 
the AF/FO order. The picture requires consistency with 



S = 1 at each V site, and with a mixture of (e£ , e£) and 
(e^aig) on each V-V pair. 

(ii) Spin wave spectra in the insulator should be consis- 
tent with exchange constants (J° J b ) set by the FO order 
(FO order is consistent with the monoclinic distortion 
involving uniform rotation of all V-V pairs in the AFI) . 
Also (see Mila et al), it can be reconciled with anomalous 
X-ray scattering results. 

2. In the PM phase 

(i) no AF order (not even a remnant of AF-LRO of 
the AFI). I-M transition strongly first order. A jump 
in the (e^e^) : ( e gi> a ig) from 2 : 1 to roughly 1 : 1 
implying a drastic rearrangement of orbital occupation 
(leading to para-orbital state?) across the I-M transition. 
The basic dependence of on orbital occupation and 
symmetry modifies these as a consequence. Recovery of 
the corundum structure in an abrupt way. 

(ii) Strong correlation driven physics as very clearly 
seen in optics and PES. Since U, V, Jh are not likely to 
vary much across the MIT, the modification of hopping in 
a way consistent with (i) holds the key to increased itin- 
erance. In any case, screening induced renormalisation of 
U, U' will occur only after the system has undergone an 
insulator-metal transition, and it is hard to understand 
how the transition itself can be "derived" by reducing 
U, U' in the PI phase. 

3. In the AFM phase 

(i) Non-FL features observed in transport studies in 
the AF-M phase, showing partial similarity to those ob- 
served in near-optimally doped cuprate superconductors 
raises interesting issues. Is this one of the elusive exam- 
ples of spin-charge separated metallic state in a three- 
dimensional oxide? 

In a multi-orbital Mott-Hubbard scenario, strong cou- 
pling to coupled orbital-spin excitations should lead to a 
dynamically fluctuating hopping, leading to inhibition of 
AF-LRO and to strongly reduced coherence, manifested 
by a low Fermi temperature. The simultaneous obser- 
vation of overdamped spin waves would also follow from 
such kind of effects. 



B. Implications for Theory 

A theoretical picture of the MIT in V2O3 must address 
these issues in a consistent way. Given that much more 
is known about the para-orbital, paramagnetic state, as 
well as the view that understanding the AFI- AFM MIT 
requires a good knowledge of the PI, we focus our at- 
tention on the para-insulating/para-metallic states. A 
complete understanding of the anomalous features near 
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the AF-I/AF-M phase is beyond current theoretical ca- 
pacity. 

Castellani et al. 8 started with a single c-axis V-V pair 
in the real crystal structure (RCS) of V2O3, and solved 
the two-site cluster including U,U' and Jh- They as- 
sumed that screening processes reduce the values of these 
parameters, and, in particular, that Jh — 0.1U. With 
this choice, and in the situation where the t 2g levels were 
split into an a\ g singlet and doublet by the trigonal 
distortion, they found that the two electrons in the a\ g 
orbitals on the pair form a total spin singlet, while the 
second electron populates the E g states. The resulting 
model is clearly a 5 = 1/2, two-orbital Hubbard model, 
with an orbital ordered, spin AF ground state. Based 
on this picture, the one-band Hubbard model was stud- 
ied extensively for twenty years with a variety of tech- 
niques. 25 As is clear from the earlier discussion, a variety 
of recent results run into direct conflict with the one-band 
modelling. 

Theoretically, the discrepancy has to do with the fact 
that J Hl which controls the spin state at each V site, is 
very poorly screened in a solid. This implies that Jh in 
V2O3 is larger than 0.1U, the value used by Castellani 
et al. Indeed, with Jh > 0.2(7, the ground state has 
been found to have 5 = 1, with a change to low-spin 
5 = state as Jh is reduced towards the value used by 
Castellani et al. 8 

Given that the occupation of the ai g ,e g orbitals 
changes discontinuously at the MIT, one would expect 
an important role for the trigonal field (since it acts like 
a fictitious external field in the orbital sector). This is 
expected to sensitively determine the occupancy of each 
orbital (orbital polarisation) in much the same way as 
the magnetisation of a paramagnet is a function of an 
applied magnetic field. In particular, one expects that 
the lower-lying orbital(s) should be more localized in the 
solid (e g ), as we shall indeed find to be the case. Further, 
the fact that the ratio of the orbital occupations changes 
discontinuously at the MIT forces one to associate a cor- 
responding change in the trigonal field as well. 

On the other hand, observation of global strong corre- 
lation signatures in various phases of V2O3 as described 
above in detail implies a fundamental inadequacy of the 
band description, and mandates use of a strong correla- 
tion picture. 

Summarising, a consistent description of the PI/PM 
MIT requires a theoretically reliable description involv- 
ing marriage of structural aspects (LDA) and strong cor- 
relation features (MO-DMFT). 

In the rest of this paper, we confine ourselves to the 
theoretical description of the PI/PM Mott transition 
in V2O3. Starting with a detailed exposition of the 
LDA+DMFT(IPT) which we use as a solver (the pros 
and cons of using IPT vis-a-vis other impurity solvers 
will be discussed), we will derive a two-fluid description 
of the PI/PM transition in V2O3 attempting to achieve 
an internally consistent description. Finally, a quanti- 
tatively accurate description of the one-particle spectral 



function across the MIT and low-T thermodynamics will 
be demonstrated within this scenario. 



III. LDA+DMFT TECHNIQUE 

As argued in detail and shown in recent work, 26 
LDA+DMFT has turned out to be the method of choice 
for a consistent theoretical description of the competi- 
tion between quasi-atomic, strong Coulomb interactions 
(multi-orbital) and itinerance (LDA spectra, encoding 
structural details in the one-electron picture) in real three 
dimensional transition-metal and rare-earth compounds. 
The central difficulty in this regard has been the choice 
of an appropriate impurity solver to solve the multi- 
orbital, asymmetric Anderson impurity problem. Two 
ways have been used with varying degrees of success: it- 
erated perturbation theory (IPT) and quantum Monte 
Carlo (QMC). 

We have used multi-orbital extension of IPT to solve 
the impurity model. On the one hand, such an approach 
should be valid if the behavior of the multi-orbital SIAM 
is analytic in U, U' , Jh'- this is known to hold for the 
general asymmetric version. MO-IPT also has the advan- 
tage of being extendable to T = 0, and the self-energies 
can be extracted at modest numerical cost. On the other 
hand, it is by no means exact, and calculations done for 
the one-band Hubbard model 26 show quantitative differ- 
ences between IPT and QMC results for the critical value 
of U = U c at which the MIT occurs. It has also been 
claimed 27 that the IPT spectral functions are very dif- 
ferent from the QMC ones, and the latter are claimed to 
be more reliable vis-a-vis the true spectral function, as 
well as with the actual, experimentally determined spec- 
tral functions. Here, we should emphasise that the IPT 
results for the many body DOS are in excellent agree- 
ment with both exact diagonalisation 28 as well as dynam- 
ical DMRG results for the one band Hubbard model in 
d = 00. While no such evidence exists for multi-orbital 
models, we believe that the above arguments show that 
IPT is a good approximation, even though it is not "nu- 
merically exact" . 

With these caveats, we describe our multi-orbital iter- 
ated perturbation theory (MO-IPT) for multi-band cor- 
related systems. For early TM oxides, one-electron band- 
structure calculations show that, in three dimensional 
cases, the t 2g DOS is well separated from the e g DOS 
as well as from the 0-2p DOS. More precisely, the "*2g" 
DOS does have contributions from components of the e g 
and 0-2p orbitals having t 2g orbital symmetry. Struc- 
tural effects, such as those produced by trigonal crystal 
fields (V2O3 29 ) and antiferroelectric distortions (VO2 30 ) 
are adequately described by LDA. In addition, the multi- 
orbital Coulomb interactions are parametrised by three 
parameters U,U' and Jh- The Hund's rule coupling, 
Jh, is very poorly screened and can be taken equal to 
its atomic value. The intra-orbital (£7) and inter-orbital 



•5 



(If) Coulomb interactions are screened in the actual 
solid: usually, their screened values have traditionally 
been calculated using constrained LDA. In correlated sys- 
tems, this is a problem, however, as the dynamical pro- 
cesses screening these parameters arise from correlated 
electrons having dualistic (itinerant-localised) character, 
rather than from free band electrons. This well-known 
problem has received scant attention to date; indeed, we 
are aware of only one previous work attempting to cure 
this malady. 31 Below, we will show how the renormalised 
U, If are self-consistently computed in a correlated ap- 
proach, and lead to a consistent description of the PES 
results in the PM phase. 

A. The many-body Hamiltonian 

Generally, the full many-body Hamiltonian for early 
TMOs is written as, 

H = ^ ( £ ka + e"(5 Qb )cj <:aCT Ck6<T + U ^ n ia]n ia l 

n ia n ib ~ Jh y ] Sj a .Sj6 (1) 

where a, b = xy,yz,zx denote the three t2 g orbitals. 
Details of the actual one-electron bandstructure in the 
real lattice structure are encoded in the one-electron 
band dispersion, e^a'. the corresponding LDA DOS is 
p(u) = N -1 Xaca ( K w _ £ ka)- Here e° a = e a - U(n aa < - 
|) + ^°~( n a<j — 1), where e a are the on-site energies of 
tig orbitals within LDA and the rest of the terms are 
subtracted therefrom in order to avoid double-counting 
of interactions already treated on the average by LDA. 
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FIG. 1. LDA partial density of states for the (red) and 
ai g (blue) orbitals, obtained from Ref. 32. 

With these qualitative remarks, we now describe our 
multi-orbital LDA+DMFT procedure. We emphasise 



that the basic method was already developed in Ref. 26 , 
and here, we extend this ideology using more detailed 
analysis to study the full one-electron spectral function in 
both insulating and metallic phases in V2O3. Our strat- 
egy is: 

(a) Beginning with LDA results in the real corundum 
lattice (see Fig. 1), derive a correlated Mott insulating 
state using multi-orbital DMFT with U = 5 eV, and 
U' = 3 eV (we use J H = 1.0 eV for U 3+ ), values ob- 
tained from constrained LDA. The LDA bandwidth is 
W = 2.5 eV, and the bare LDA trigonal field is read off 
as A = 0.32 eV. In what follows, we will work in the ba- 
sis of LDA eigenstates which diagonalises the one-particle 
density matrix. 

(b) Mimic the effects of external pressure by noticing 
that it should lead to modification of the renormalised 
(correlated) value of the trigonal field. In line with this 
ideology, first known to be propounded by Mott and co- 
workers, we search for the instability of the correlated 
(Mott insulator) solution found in (a), to a second so- 
lution of the DMFT equations as a function of A. We 
emphasise here that we do not change the bare LDA pa- 
rameters: indeed, we argue that a Mott transition from a 
correlated insulator to correlated metal cannot be validly 
described by changing bare LDA parameters, since these 
have no clear meaning in a strongly correlated system. 

(c) To provide a quantitative description of the one- 
electron spectral function in the metallic phase, we use 
the correlated (DMFT) results to compute the screening- 
induced reduction in U, U' in the metallic phase. This is 
crucial: we derive the screened U, If in the PM phase 
after deriving the I-M transition, and do not derive the 
I-M transition itself by reducing U, If, as seemingly done 
in earlier work. 26,32 Using the DMFT result, the screened 
U, If are estimated by an extension of Kanamori's t- 
matrix calculation to finite density. 

Using the screened values of U, If (notice that Jh is 
almost unaffected by screening, so we use the same value 
for it throughout), we compare our theoretical (corre- 
lated) DOS with PES and XAS results obtained experi- 
mentally in the PM phase. 

Incorporation of electron correlations into the LDA 
gives rise to a two-stage renormalisation: 

(1) U, If and Jh give rise to multi-orbital Hartree 
shifts in the on-site orbital energies of each ti g orbital. 
In V2O3, the trigonal field splits the ti g degeneracy, with 
the lowest a\ g orbital [~ (xy + yz + zx)\ lying about 
A = 0.32 eV below the higher lying orbitals [~ (xy — 
yz), (2xy — yz — zx)\ within LDA. Given this, the a\ g 
orbital is always occupied by one electron, the second 
residing in the e£ orbitals. The observation of S = 1 on 
each V site requires strong Jh, implying even stronger 
U, U' , even in the PM phase. 

Multi-orbital Hartree shifts renormalise the orbital en- 
ergies: e aig = e + U'n e * (e e j = e + A + U'n aig ), where 
n e j (n 0lB ) is the (dig) orbital occupation. These 
shifts correspond to effects captured by LDA+U. 12 They 
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do give the correct, insulating ground states (with or- 
bital/magnetic order), but cannot describe the phase 
transition(s) from correlated Mott insulators to corre- 
lated metals. This can be traced back to the fact that 
LDA+U treats correlations on a static level, neglecting 
quantum nature of electron dynamics, and so cannot ac- 
cess the spectral weight transfer-driven physics at the 
heart of Mott-Hubbard transitions. 

(2) In a one-electron picture, this would be the end of 
the story. In reality, however, hopping of an electron from 
a given site to its neighbor(s) is accompanied by dynami- 
cal generation of particle-hole pairs (the more, the larger 
U, U' are), which inhibit its free band motion. Electrons 
can move quasicoherently by dragging their correspond- 
ing "electronic polarisation cloud" along. With increas- 
ing U, V, electrons get more and more "localised", cor- 
responding to transfer of coherent low energy spectral 
weight to high-energy (quasi-atomic) incoherent regions, 
until at the MIT, all the weight resides in the incoherent 
Mott-Hubbard bands. It is precisely this effect that is 
out of scope of LDA+U, and requires dynamical mean 
field theory (DMFT) for a consistent resolution. 

Since the system is strongly correlated, the small 
changes in bare LDA parameters caused by (1) lead to 
large changes in transfer of dynamical spectral weight. 
Specifically, in systems undergoing MIT, small changes 
in bare LDA values of lattice distortion(s) transfer high- 
energy spectral weight to low energies, driving the Mott 
transition. 



B. The one-particle Green's functions 

Given the actual LDA DOS for the t 2g orbitals (this 
includes the V-d orbitals and O 2p part having u t 2g " 
symmetry), the band Green's functions within the LDA 
(in the basis which diagonalises the one-particle density 
matrix) are G ab (uj) = 5 ab G a {uj) = S a bjj Ek( w _ £ka) _1 - 
We define the correlated one-electron Green's function 
and the associated irreducible self-energy for each orbital 
a, by G acr (cL>) and £ OCT (aj): the two are related by the 
usual Dyson's equation, 



(2) 



It is obvious that the Green's functions can be exactly 
written down for the non-interacting case, as well as for 
the atomic limit (eka = 0) . In contrast to the case of the 
one-band Hubbard model, however, the exactly soluble 
atomic limit contains the local, inter-orbital correlation 
function, < n a nf, >, in addition to < n a >. 



1. MO-IPT: an interpolative ansatz for multi-orbital systems 

In the spirit of the IPT developed by Rosenberg et 
al. for the one-orbital Hubbard case, we require an in- 
terpolative scheme that connects the two exactly soluble 



cases above, gives correlated Fermi liquid behavior in the 
metallic phase, and a Mott-Hubbard transition from a 
correlated FL metal to a Mott insulator as a function of 
U, U' for commensurate cases. 33 In order to achieve this, 
we have extended the philosophy of Rcf. 35. The cen- 
tral requirements for a consistent interpolative scheme 
capable of describing all of the above are that: 

(i) Formally defined one-electron Green's function, 



(3) 



where e^ a describes the dispersion of the LDA bands for 
orbitals a,b = t 2g , and the self-energy is given by 



S„(w) 



(4) 



with 



E^M = N a & J2 Gl^Gliiw^GHiwt + uo m - zlu) 

' lm 

(5) 

being the second-order (in U, U') contribution. Here, 
N a b = 2 for a, b = e*i,e* 2 an d 4 for a,b = ai g ,e g i t2 - 
Finally, the bath propagator is given as 



G» = 



UJ + H a - A a (uj) 



(6) 



with A a (w) interpreted as the dynamical Weiss field for 
orbital a. 

(ii) The interpolative self-energy for each orbital a 
should be chosen by fixing interpolative parameters such 
that the exact Friedel-Luttinger sum rule is strictly (nu- 
merically) obeyed, and, 

(Hi) to reproduce the Mott insulator beyond a criti- 
cal coupling, a high-energy expansion around the atomic 
limit is performed, yielding another equation for the in- 
tcrploative parameters. Here, the high-energy expan- 
sion is truncated by including only the first few terms 
which guarantee the exact reproduction of the first three 
moments of the one-electron spectral function. In con- 
trast to the one-band case, however, the exact atomic 
limit for the multi-orbital case contains the local, inter- 
orbital correlation function, D ao \n\ =< n a ni, >. 34 We 
are aware of only one earlier work 35 where -D a f,[n] is com- 
puted using the coherent potential approximation (CPA). 
Strictly speaking, this is an approximation to the Hub- 
bard model(s) which is qualitatively valid in the Mott 
insulating state, but is known to fail in the correlated 
PM phase(s). This is because CPA replaces the actual, 
dynamical (annealed) "disorder" in the PM phase(s) by 
quenched, static disorder, and thus fails to capture the 
dynamical Kondo screening central to deriving correct 
(correlated) FL behaviour in the PM phase. Given this, 
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it is hard to identify the extent to which computed re- 
sults depend upon introducing such approximations, and 
this should be checked carefully by comparison with cal- 
culations which compute all local correlators in a single, 
consistent scheme. The correct way to compute D a b[n] is 
actually not complicated within multi-orbital IPT, and 
is described below. 

These two equations for the parameters A a t> and B a b 
are solved to yield these as explicit functions of U,U',< 
n a >,< n a n,b > and < n° > (this last average is the "ef- 
fective" number of fermions in orbital a corresponding to 
an "effective" Green's function used in the interpolative 
IPT, see Refs. 35, 36). Explicitly, we have, 



read off directly from Fig. 2. The orbital occupations are 
computed to be {n aig , n e * , n e ^ — 0.32,0.34.0.34) in the 
PI, in nice agreement with XAS estimations. 



and 



A ab = 



B„ 



n a (l - 2n a ) + Dgb[n\ 

»g(i-»2) 

(1 - 2n a )U ab + fx - n a 
C>2(l-n2) 



(7) 



(8) 



where n a and rP a are defined from the GFs G a (oS) and 
G°(lo). The inter-orbital correlation function D a b[n] is 
calculated from 



1 r°° 

Dab[n] =< n a >< n b > +— / Z(w)f{uj)du) , (9) 

U ab J ~ oo 

the last term following from the equation of motion for 
G«(w) and S(w) = -±Im [E (w)G„(w)]. 

The above equations form a closed set of coupled, 
non-linear equations which are solved numerically. We 
found fast convergence of the self-consistent system of 
equations, and typically twenty iterations sufficed for the 
parmeter values considered here. The converged results 
allow us to study the one-particle DOS, and the corre- 
sponding orbital occupations, spin states, as well as the 
strength and character of local multi-orbital correlations 
in both PI and PM phases, as described below in detail. 



IV. RESULTS AND DISCUSSION 

In this section, we present a detailed set of results for 
the one-particle spectral function in both the PI and PM 
phases in V2O3. While doing so, we will make extensive 
contact and discuss important differences between our 
work here and previous results recently obtained by other 
authors. 26,32 

In Fig. 2, we show the single particle DOS for the Mott 
insulator, obtained with U = 5.0 eV, U' = 3.0 eV and 
Jh = 1.0 eV as correlation parameters for this system. 
These are slightly different from those used in our previ- 
ous work, 29 but are roughly the same as those used by 
Held et al. recently. 32 A clear Mott-Hubbard gap, Eq — 
0.2 eV, is seen, and, as expected from the orbital assign- 
ment, the states are more localised in the solid. The 
renormalised trigonal field A[ — S aig — 5 e -* = 0.32 eV, is 




FIG. 2. Orbital-resolved (upper panels) and total (lower 
panel) one electron spectral function for the insulating phase 
of V2O3 obtained with 5.0 eV and Jh = 1.0 eV. 

We now study the paramagnetic metallic state ob- 
tained as an instability of the correlated Mott insulator 
derived above. In other works, the PM state is "derived" 
not by searching for an instability of the correlated PI 
state under pressure, but by computing the LDA band- 
structure for the "metallic" state without correlations. 
The screened values of U, U' are computed using con- 
strained LDA, and these are then used to describe the 
PM phase. In reality, however, one has to study the 
transition to the PM phase without leaving the corre- 
lated picture, and derive the transition by searching for 
the second, metallic solution of the DMFT equations un- 
der pressure. 

To justify our new approach, we specify the problems 
associated with earlier approaches: 26 ' 32 ' 37 

(1) It is theoretically inconsistent to derive a phase 
transition between two strongly correlated phases by us- 
ing corresponding LDA bandstructures to separately de- 
rive the two phases. This is because using changes in bare 
LDA parameters to study correlated phases is clearly 
problematic, since these parameters have no clear mean- 
ing in a correlated picture. One must use the renor- 
malised values of these parameters instead, and these are 
generically modified in unknown ways by strong multi- 
orbital correlations. These changes in bare LDA param- 
eters, and the modification of the response of correlated 
electrons to these changes, must be selfconsistently de- 
rived within the LDA+DMFT procedure. Clearly, this 
route has not been used in other approaches. 

(2) It follows that an inescapable consequence of us- 
ing such approaches is that the values of U, U' used for 
the PM phase are computed using constrained LDA (i.e, 
using the uncorrelated bandstructure assuming that the 
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screening electrons are free band electrons). However, 
in reality, the screening electrons in the correlated PM 
phase have a dualistic character generic to the Mott- 
Hubbard character of the system. As is known from 
Ref. 25 , the electronic kinetic energy, or itinerance, is re- 
duced in the PM phase: it is these correlated electrons 
which screen U, V in the real correlated system. In 
an "ab initio" treatment, the effective U, U' should be 
computed using the correlated spectral functions to esti- 
mate screening. Replacement of the renormalised spec- 
tral functions by bare LDA ones will introduce an ap- 
proximation, overestimating the screening of U, V (this 
is hard to quantify, but is estimated to be of order of 
twenty percent!). 

In order to avoid these difficulties, we adopt the fol- 
lowing strategy. 

(A) We hypothesise that external pressure modifies the 
trigonal field. To our knowledge, this is not completely 
new: Mott and co-workers proposed such ideas in the 
seventies, 1 and more recently, Tanaka made a similar hy- 
pothesis in a cluster approach for V2O3. 11 To model this 
change in A t under pressure, we do not change the trig- 
onal field by hand. Rather, we input trial values of A t , 
changing it from its value in the (Mott) PI by small trial 
amounts, and search numerically for its critical value, A£ , 
which stabilises the second, correlated metallic solution 
of the DMFT equations. The new values of At in the 
(correlated) PM phase are again read off from the con- 
verged DOS for each orbital. We emphasise that we do 
not decrease U, U' by hand, neither do we use different 
LDA DOS for different phases, for reasons explained be- 
fore. We note that Savrasov et al. 38 have employed simi- 
lar ideology to study the giant volume collapse across the 
a — 5 transition in Pu. 

In a strongly correlated system, small changes in the 
renormalised trigonal field lead to large changes in dy- 
namical spectral weight transfer from high- to low en- 
ergies, typically over a scale of a few eV. This is pre- 
cisely our mechanism for the first-order Mott transition 
in V2O3 under pressure. We expect the free energy to 
have a double well structure. Pressure changes the trig- 
onal field (we remind the reader that A t acts like an 
external field in the orbital sector), lowering the second 
minimum (PM) relative to the first (PI) beyond A£. 

(B) Using the converged DOS for each orbital, the oc- 
cupation^) of various orbitals (and their changes from 
their PI values), the local spin value at each V site, and 
information about the detailed character of the PM state 
is directly obtained. 

In Fig. 3, we show our results for the PM phase ob- 
tained within our technique. At T = 0, the hypothet- 
ical PM phase (it is never observed in reality) shows a 
sharp, quasicoherent FL resonance. We identify this fea- 
ture with combined spin-orbital Kondo screening in the 
PM phase of a multi-orbital Hubbard model. This is eas- 
ily seen as follows. To obtain a correlated Mott insulator, 
we need not only U = 5.0 eV, but also U' ~ (U — 2Jh) = 
3.0 eV: indeed, if U' were ignored, a t2 g electron hopping 



from one V site to its neighbor could always hop off like 
a band electron just by going into an unoccupied t2 g or- 
bital at that site, making a PI state impossible. Given 
the small spectral weight carried by this feature, we ex- 
pect a low lattice coherence scale, above which the PM 
would be described as an incoherent metal. The trigonal 
field in the PM, A' t = —0.291 eV and the occupations 
of each orbitals, [n i s ,n es i,n efl 2 = 0.38,0.31,0.31], 
are read off from the converged PM solution of the DMF 
equations. Very satisfyingly, the spin state remains un- 
changed, and the orbital occupations change across the 
MIT in semiquantitative agreement with XAS results. 10 
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FIG. 3. (color online) Orbital-resolved (upper panels) and 
total (lower panel) one electron spectral functions for the 
metallic phase of V2O3. Note that only the ai g orbital DOS 
crosses Ef in the metallic phase; the orbitals still shows 
Mott-Hubbard insulating features, showing the "two-fluid" 
character of the MIT in V2O3. 

In Figs. 4 and 5, we show the effect of finite temper- 
ature on our results. As expected on general grounds 
within the DMFT framework, the FL resonance is broad- 
ened by finite-T and lowered in height (the pinning of 
the interacting DOS at Ep to its LDA value, dictated by 
Luttinger's theorem, holds only at T = 0). The effects 
of introducing chemical disorder in the PI is shown in 
Fig. 6. The results were obtained by combining multi- 
orbital IPT with the coherent-potential approximation 
(CPA) 39 . In agreement with very recent observations, 23 
we indeed observe a broadened "quasiparticle" in the 
PM, and closing in of the Mott gap in the PI by in- 
coherent spectral weight transferred across large energy 
scales from high- to low energies. Comparing Figs. 4 
and 6, it is clear that, at sufficiently high-T, the spectra 
in the chemically disordered PI and the PM phases do re- 
semble each other qualitatively. As observed by Allen, 40 
this implies that there is no fundamental difference be- 
tween the "metal" and "insulator" at sufficiently high T: 
this agrees with the observation that the first order Mott 
transition is replaced by a smooth crossover at high T . 
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FIG. 4. Effect of temperature (T) on the orbital-resolved 
(upper panel) and total (lower panel) one electron spectral 
functions for the insulating phase of V2O3. 





FIG. 5. Effect of temperature (T) on the orbital-resolved 
(upper panels) and total (lower panel) one electron spectral 
functions for the metallic phase of V2O3. 

An extremely important conclusion follows directly 
from an examination of the orbital-resolved DOS in the 
PM phase. We find that the e£ orbital DOS shows "Mott 
insulating" (see Fig. 3) behavior, while only the a\ g or- 
bital DOS is responsible for the metallicity. This con- 
stitutes an explicit realisation of the "two-fluid" model 
used phenomenologically in connection with the MIT in 
disordered semiconductorsin the past. 1 In Refs. 29,30 we 
already showed the orbital selective character, as well 
as the evolution of the DOS at Ep as a function of the 
occupation of the a\ g orbital. A clear first-order I-M 
transition around n ai = 0.38 was found, involving, as 
described above, a discontinuous change in (selfconsis- 
tently determined) occupations of each orbital. These 
observations are intimately linked to the multi-orbital 



Mott-Hubbard character of correlations in V2O3. Po- 
larised XAS results might already hold the clue to es- 
tablishing an approximate two-fluid character of the PM 
phase: the a\ g spectral weight should dominate over the 
contribution for energies up to the Mott gap. Orbital 
resolved optical studies could also be used to test our 
picture. 




FIG. 6. Effect of disorder (v) on the orbital-resolved (upper 
panels) and total (lower panel) one electron spectral functions 
for the metallic phase of V2O3. 

Strong indirect support for our picture comes from the 
early observation 9 of an anisotropic change in the lattice 
constants along a/b (planar) and c axes across the P-MIT 
in V2O3. Instead of a uniform volume collapse expected 
across the MIT, 41 increase in a(b) and a decrease in c was 
found across the MIT. Such an anisotropic volume change 
across the MIT is inconsistent with simultaneous gapping 
of all t2g orbitals (where we would expect an isotropic 
volume change), but is completely consistent with our 
(orbital selective) two-fluid picture derived above. 



V. COMPARISON WITH PES AND XAS 

In this section, we describe how our approach provides 
an excellent description of the experimental photoemis- 
sion (PES) and X-ray absorption (XAS) data on V2O3 in 
the PM phase. As argued before, this requires us to re- 
compute the full one-particle local spectral function (to- 
tal DOS) using values of U, U' renormalised by dynamical 
metallic screening in the correlated metallic phase. In or- 
der to do this, we have used an extension of Kanamori's 
i-matrix approach 3 to estimate these parameters. 

In the multi-orbital case, this is a horrendous problem 
in general. Fortunately, in the effective two-fluid picture 
of the PM phase derived above, the general analysis can 
be simplified. This is because the e g electrons remain 
"insulating" , i.e, Mott localised, up to energies of the or- 
der of the Mott-Hubbard gap. We then expect only the 
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a,\g electrons to provide efficient screening, and so con- 
sider only the a\ g band in the computation of the effective 
U, U' below. In general, we need the full q-dependent 
particle-particle susceptibility for this purpose. Using the 
LDA+DMFT Green function for the a lfl orbital, 

Xpp(q) = -J7 ^2 Ga ls (q-k,iw m -iv n )G aig (\L,iv n ) . 

nmk 

(10) 

In d = oo, this can be expressed as an integral over 
the LDA DOS and the full irreducible one-electron self- 
energy 25 , permitting a direct evaluation. The onsite 
Hubbard U is renormalised by the local part of this sus- 
ceptibility, via the equation, 



U, 



u 



eff 



i + u Xlo M 



(11) 



Using the relation U ~ (U' + 2Jh), valid for t2 g 
systems, along with the fact that Jjj is essentially un- 
screened, we estimate U, U' in the PM phase. We observe 
that this implies a frequency-dependent U e ff = U{ui). 
We have found, however, on computation that the in- 
dependence is weak for energies up to the Mott gap, and 
so use its uj = value U e ff = U(0) in what follows. 
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FIG. 7. (color online) Comparison of theoretical LDA 
+DMFT result (blue) for the total one-electron spectral func- 
tion in the metallic phase of V2O3 to the experimental results 
taken from Refs. 21, 22 (for PES) and from Ref. 42 (for XAS). 



clearly observed. In addition to the detailed shape of 
the lower Hubbard band (in PES), excellent agreement 
with the intense peak in XAS is also clear. Consider- 
ation of parts of the spectrum for lu < —3.0 eV and 
lu > 1.2 eV is hampered by our restriction to the ti g 
sector in the LDA+DMFT calculations. Due to the re- 
duction of U, V as above, the ti g orbital occupation is 
now (n a i g ,n eg i,n eg 2) = (0.36,0.32,0.32), in even better 
agreement with XAS results. 10 

However, though good, the agreement is not quite so 
perfect in the low-energy region: our computed '"broad" 
peak (ascribed to a "quasiparticle" in earlier work 22 ) is 
narrower than the experimental feature by a factor of 
f .8. On first sight this might seem to confirm the inter- 
pretation in the earlier work. However, we observe that 
this feature is peaked at u> = —0.37 eV, while a clear 
pseudogap-like dip is resolved around Ep. Hence, in our 
picture, the metallic phase cannot be described in a FL 
quasiparticle language; instead, short-lived, incoherent, 
non-FL pseudoparticles should dominate the PM phase. 
Interestingly, observation of a linear-in-T(instead of the 
T 2 form for a correlated FL) resistivity supports a non- 
FL quasiparticle interpetation. It is possible that the T 
regime where this is valid lies above an effective FL co- 
herence scale (below which a T 2 term in resistivity would 
follow) which is masked by emergence of orbital/spin or- 
dered insulating states at lower T. At T > T co h, the 
dc resistivity is indeed linear in T in a Hubbard model 
framework, where it arises from inelastic scattering off 
unquenched spin-orbital local moments in a d — 00 multi- 
band Hubbard model. 

Our observation of a low energy pseudogap feature can 
be traced back to the strong bonding-antibonding split- 
ting observed in LDA results (see Fig. f). This is a di- 
rect consequence of strong hopping along the a\ g orbitals, 
leading to strong covalency and robust singlet character 
between V-V pairs along the c-axis. (see Ref. 43 for an 
early discussion on this point). Our analysis does par- 
tially show up the effects of strong covalency, manifested 
in the pseudogap feature found above. However, the fact 
that we can resolve most of the spectrum accurately, but 
fail to reproduce the correct broadening of the low-energy 
feature, implies that it may be necessary to explicitly 
consider the dynamical effects of intersite (V-V) corre- 
lations for a complete resolution of the PES spectrum, 
as alluded to in the experimental section. This is how- 
ever out of scope of LDA+DMFT, and requires a cluster 
extension. 



Starting with the values of U, U' used earlier, we es- 
timate x'ioc(°) - °- 084 > yielding U e}s ~ 3.5 eV. With 
Jh = f.O eV, this implies that U' ~ f.5 eV. We have 
recomputed the one-electron spectral function for the 
PM phase using these values. The results are com- 
pared with experimental work 21 ' 22 ' 42 in Fig 7. Very sat- 
isfyingly, excellent quantitative agreement over almost 
the whole energy scale from —3.0 < co < 1.2 eV is 



VI. CONCLUSION 

In conclusion, we have studied the first-order Mott 
transition under pressure in V2O3 using the state-of-the- 
art LDA+DMFT technique. We have proposed a new 
picture for the MIT, which is driven by large changes in 
the transfer of dynamical spectral weight (via DMFT) 
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accompanying small changes in the renormalised trigo- 
nal field splitting under pressure. Very good quantitative 
agreement with the orbital occupations, spin state of V 3+ 
ions, as well as effective mass enhancement in the PM 
state is obtained. The MIT is found to be first-order, and 
orbital selective (only the a\ g DOS shows metallic behav- 
ior). Finally, using the correlated solution, we have com- 
puted the screening induced renormalisation of U, U' in 
the PM phase. Using these, excellent quantitative agree- 
ment with the full one-particle spectral function (PES 
and XAS) is found in the PM phase. These findings con- 
stitute strong support for our underlying two-fluid pic- 
ture, which is ultimately an interesting manifestation of 
strong, multi-orbital Coulomb interactions in this early 
transition-metal oxide. 
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